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Abstract 
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other classes of functions. 
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1 Introduction 



Pseudo-differential equations for complex-valued functions defined on a non- Archimedean lo- 
cal field are among the central objects of contemporary harmonic analysis and mathematical 
physics; see the monog raphs [Vladimirov et al. 1994] . [Kochubei 2001] . (Alb everio et al. 2010] , 
and the survey [Zuniga-Galindo 2011 . 

The simplest example is the fractional differentiation operator a > 0, on the field Qp of 
p-adic numbers (here p is a prime number). It can be defined as a pseudo-differential operator 
with the symbol |^|p where | ■ |p is the p-adic absolute value or, equivalently, as an appropriate 
convolution operator. 

Already in this it was first shown by Vladimirov (see [Vladimirov et al. 1994] ). 

properties of the p-adic pseudo-differential operator are much more complicated than those 
of its classical counterpart. It suffices to say that, as an operator on L2(Qp), it has a point 
spectrum of infinite multiplicity. Considering a simple "formal" evolution equation with the 
operator D"' in the p-adic time variable t, Vladimirov (see [Vladimirov 2003] ) noticed that such 
an equation does not possess a fundamental solution. 

At the same time, it was found in [Kochubei 2008] that some of the evolution equations 
of the above kind behave reasonably, if one considers only solutions depending on \t\p. This 
observation has led to the concept of a non-Archimedean wave equation possessing various 
properties resembling those of classical hyperbolic equations, up to the Huygens principle. 

In this paper we consider the Cauchy problem for a class of equations like 

D''u + a{\x\p)u = f{\x\p), a;eQp, (1.1) 

assuming that a solution is looked for in the class of radial functions, u = u{\x\p)] the accurate 
definition of and assumptions on a, / are given below. This Cauchy problem is reduced to 
an integral equation resembling classical Volterra equations. It appears that the equation (1.1) 
and its generalizations considered on radial functions constitute p-adic counterparts of ordinary 
differential equations. 

The author is grateful to M. M. Malamud for a helpful discussion. 

2 Preliminaries 

2.1. Local fields. Let if be a non-Archimedean local field, that is a non-discrete totally 
disconnected locally compact topological field. It is well known that K is isomorphic either to 
a finite extension of the field Qp of p-adic numbers (if K has characteristic 0), or to the field of 
formal Laurent series with coefficients from a finite field, if K has a positive characteristic. For 
a summary of main notions and results regarding local fields see, for example, [Kochubei 2001] . 

Any local field K is endowed with an absolute value \ ■ \k) such that \x\k = if and only 
if a; = 0, \xy\K = \x\k ■ \y\K, \x + y\K < max(|a;|ii-, \y\K)- Denote O = {x & K : \x\k < 1}, 
P = {x & K : \x\k < Oisa subring of and P is an ideal in O containing such 
an element (3 that P = (30. The quotient ring 0/P is actually a finite field; denote by q its 
cardinality. We will always assume that the absolute value is normalized, that is \P\k = 
The normalized absolute value takes the values , N ^ Z. Note that for if = Qp we have 
(3 = p and q = p', the p-adic absolute value is normalized. 
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Denote by S* C O a complete system of representatives of the residue classes from 0/P. 
Any nonzero element x & K admits the canonical representation in the form of the convergent 
series 

X = (xo + + 0:2/3' + ■■■ ) (2.1) 

where n G Z, \x\k = q"", Xj E S, Xq ^ P. For K = Qp, one may take S = {0,1, . . . ,p — 1}. 

The additive group of any local field is self-dual, that is if x is a fixed non-constant complex- 
valued additive character of K, then any other additive character can be written as Xa{x) = 
x{(ix), X E K, for some a & K. Below we assume that x is a rank zero character, that is 
x(x) = 1 for X G O, while there exists such an element xq E K that \xo\k = q and x{.xq) 7^ 1. 

The above duality is used in the definition of the Fourier transform over K. Denoting by 
dx the Haar measure on the additive group of K (normalized in such a way that the measure 
of O equals 1) we write 

m = j xixOf{x)dx, ieK, 

K 

where / is a complex- valued function from Li{K). As usual, the Fourier transform J-" can be 
extended from Li{K) fl L2{K) to a unitary operator on L2{K). If J-"/ = / G Li{K), we have 
the inversion formula 

f{x) = J xi-xomd^- 

K 



2.2. Integration formulas. As in the real analysis, there are many well-known formulas 
for integrals of complex-valued functions defined on subsets of a local field. There exist even 
tables of such integrals [Vladimirov 2003] . Note that formulas for integrals on Qp and its 
subsets, as a rule, carry over to the general case, if one substitutes the normalized absolute 
value for I ■ Ip and q for p. 



Here we collect some formulas used in this work. 



1 — q~^ 

\x\'^'^ dx = —q""", here and below n E Z,a > 0. 



\x\jf<q" 



1-q- 



I X Qi I dx 



g-2 + g-", 
g(l - q- 



\0'\k-> — q 



\x\K=q" 



log \x\k dx 



n 



q-l 



|a;k<g" 
log \x — a\K dx 



log \a\ 



K 



(T logg. 

log? 
q-\ 



\(i'\K, \(i'\K — q 



\x\K=q" 



dx = g"; 



dx 



1 - - ) g"- 

1/ 



\x\K=q" 



(2.2) 
(2.3) 
(2.4) 
(2.5) 
(2.6) 
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dx = q''-\ Oy^koeS (2.7) 
\^\K=q",xo=ko 

(the restriction xq = ko is in the sense of the canonical representation (2.1)). 

dx= (^l- 1^ g". (2.8) 

\x\K=q", ^oT^ko 

2.3. Test functions and distributions. A function / : — C is said to be locally 
constant, if there exists such an integer / that for any x & K 

f{x + x') = /(x), whenever \x'\ < q~K 

The smallest number / with this property is called the exponent of local constancy of the 
function /. 

Let T>{K) be the set of all locally constant functions with compact supports; it is a vector 
space over C with the topology of double inductive limit 



V{K) = Ih^ V^^ 



where Pjy is the finite-dimensional space of functions supported in the ball = \^x ^ K : 
\x\ < q^} and having the exponents of local constancy < I. The strong conjugate space T>'{K) 
is called the space of Bruhat-Schwartz distributions. 

The Fourier transform preserves the space V{K). Therefore the Fourier transform of a 
distribution defined by duality acts continuously on V{K). As in the case of R", there exists 
a well-developed theory of distributions over local fields; it includes such topics as convolution, 
direct product, homogeneous distributions etc (see [Vladimirov et al. 1994] . [Kochubei 2001] . 
[Albeverio et al. 2010] ). In connection with homogeneous distributions, it is useful to introduce 
the subspaces of V{K): 

<ij(K) = G V{K) : V^(0) = 0}, 

$(ir) = jv? G P(ir) : J^{x)dx = 

The Fourier transform J-" is a linear isomorphism from \E'(i^) onto ^{K), thus also from 
^'{K) onto "^'{K). The spaces $(if) and ^'{K) are called the Lizorkin spaces (of the second 
kind) of test functions and distributions respectively; see jAlbeverio^ et al. 2010] . Note that two 
distributions differing by a constant summand coincide as elements of ^'{K). 

3 Fractional Differentiation and Integration Operators 

3.1. Riesz kernels and operators generated by them. On a test function ip G V{K), the 
fractional differentiation operator D°', a; > 0, is defined as 

(DV)(x)=.F-M|eiV(^))(e)](x). (3.1) 
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However does not act on the space V{K), since the function ^ i— )■ |^|" is not locally constant. 
On the other hand, : (^{K) ^ ^{K) and : ^'{K) ^'(K); see [Albeverio et al. 2010] . 
and that was a motivation to introduce these spaces. 

The operator can also be represented as a hypersingular integral operator: 

(^V) (x) = Y^^^ / \yr~'Hx -y)- ^(x)] dy (3.2) 

K 

[Vladimirov et al. 19941 IKochubei 200T] . In contrast to (3.1), the expression in the right of 
(3.2) makes sense for wider classes of functions. Below we study this in detail for the case of 
radial functions. 

27rz 

The expression in (3.2) is in fact a convolution /_q,*v^ where /s, s G C, s ^ 1 (mod Z), 

logg 

is the Riesz kernel defined first for Re s > 0, s 7^ 1, as 

^ ' Tk^sY ^ ' 1 - g-^ ' 
and then extended meromorphically, as a distribution from V{K) given by 



1 — q ^ 1 — q 



\^\k -J \'^\k 

x\k>1 |a;|if<l 



27r^ 27rz 

s 7^ 0, s ^ 1 (mod Z). For s = 0, /o(x) = 5{x). For s = 1 (mod Z), we define 

log q log q 



fs{x) = ^log|x| 

logg 



It is well known that fs* ft = fs+t in the sense of distributions from V{K), ifs,t,s + t^l 
2Tii 

(mod Z). If these kernels are considered as distributions from ^'(K), then fs*ft = fs+t for 

logg 

all s,t eC [Alb everio et al. 2010] . In view of this identity, it is natural to define the operator 
D"", a > 0, setting 

(x) = (/« * v^)(x) = j\x- yrxMy) dy, V e V{K), a ^ 1, (3.3) 

K 

and 



1 — g f 

(z^^V) {x) = -j — / log k - y\Kv{y) dy. (3.4) 
g log g y 



K 



Then D°'D^° = / on V{K), if a 7^ 1. This property remains valid on ^{K) also for a = 1. 

The above notions and results are well known; see [Vladimirov et al. 1994j . 
[Albeverio et al. 2010] . We come to new phenomena considering the case of radial functions. 



3.2. Operators on radial functions. Let m be a radial function, that is m = M(|x|i^ 

l^q 

1-q 

X G K, we denote by Xq the element from S C O appearing in the representation (2.1 



X G K. Let us find an explicit expression of W^u, a > 0. Below we write da = -a-i ' "^^^ 
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Lemma 1. If a function u = u{\x\k) is such that 

m oo 

J2 < oo, \u{q')\ < oo, (3.5) 



fe=— oo l=m 



for some m G Z, then for each n G Z the expression in the right-hand side of (3.2) with 
(p{x) = u{\x\k) exists for \x\k = q"', depends only on \x\k, and 



1 \ °° 

- J] q--^u{q^). (3.6) 
;=n+i 



fe = — CXD 

+ da ( 1 

Proof. We find, using the ultrametric properties of the absolute value, that 
iD-u){x) = J [ui\x - vIk) - u{\x\k)] dy. 

\y\K>\x\K 

If \y\K = \x\k and yo ^ Xq, then the integrand vanishes. Therefore by (2.6), 

n-l 



/fc J- n 

(D'^«)(x) = ^ y \x\-^"-' [u{q') - u{q^)] dy 

k=—oo, I t 

oo „ 

^djl- -] ^ g'^ [u{q'') - u{q^)] +djl--) g""' [^(gO - «(g")] ■ 

^ fe=-oo ^ l=n+l 

It is clear from this expression that {D°'u){x), \x\k — g", depends only on \x\k- After 
elementary transformations we get (3.6). ■ 

Definition. We say that the action D^u, a> on a radial function u is defined in the strong 
sense, if the function u satisfies (3.5), so that D°'u{\x\k), \x\k 0, is given by (3.6), and there 
exists the limit 

D^uiO) '^^^ ]im D'^u{\x\k). 

It is evident from (3.2) that annihilate constant functions (recall that in ^'{K) they are 
equivalent to zero). Therefore D~°' is not the only possible choice of the right inverse to D"'. 
In particular, we will use 

(/»(x) = (i^-»(x) - (i^-»(0). (3.7) 
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This is defined initially for (p e V{K). It is seen from (3.3), (3.4), and the ultrametric property 
of the absolute value that 

(^»(^) = ^f^ / (k-y|r'-|y|r')(^(y)dy, a ^ 1, (3.8) 
\y\K<\AK 

and 

1 — q f 

{I^ip){x) = — / {\og\x -y\K -log \y\K)(p{y)dy. (3.9) 

glogg ,/ 

\y\K<\'MK 

In contrast to (3.3) and (3.4), in (3.8) and (3.9) the integrals are taken, for each fixed x E K, 
over bounded sets. 

Let us calculate I°'u for a radial function u — u{\x\k)- Obviously, {l'^u){0) — whenever 
7° is defined. 

Lemma 2. Suppose that 

m 

J2 max g"'^) |m(/) I < oo, if a ^ 1, 

k=—oo 

and 

m 

J2 \k\q'' \u{q'')\ < oo, if a = 1, 

for some m G Z. Then I^u exists, it is a radial function, and for any x ^ 
{ru){\x\K) = q-''\xrM\x\K) + V"i / {\^rK'-\y\K~')^{\y\K)dy, a^l, (3.10) 



and 



1 - g° 

\y\K<\x\K 



1 — q f 

{I^u){\x\k) = q~^\x\Ku{\x\K) + —, / {log\x\K -log \y\K)u{\y\K) dy. (3.11) 

g log g J 

\v\k<\x\k 

Proof It is sufficient to compute the integrals over the set {y E K : \y\K — \x\k}, and that 
is done using the integration formulas (2.3) and (2.5). ■ 

It follows from Lemma 2 that the function I°'u is continuous if, for example, u is bounded 
near the origin (see an estimate of the integral /„.o in the proof of Theorem 1 below). If 
li^dajlii:)! < C|x|^^, as \x\k > 1, then \ {I°'u) {\x\k)\ < C'|a;|^~^, as \x\k > 1- Here and below we 
denote by C various (possibly different) positive constants. 

It is easy to transform (3.10) and (3.11) further obtaining series involving u{q^). 

Obviously, D"7" = 7 on V{K), if a 7^ 1, and on $(7^), if a = 1. Since by Lemma 1 and 
Lemma 2, the operators are defined in a straightforward sense for wider classes of functions, it 
is natural to look for conditions sufficient for this identity. 
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Lemma 3. Suppose that for some m G Z, 

m oo 

J2 max(g^g"'=) < cx), J] < oo, 

fc=— oo l=m 

if a 1, and 

m oo 

J] Ifclg'^ < OO, ^ / |i;(g') | < oo, 

fc=— CO l=m 

if a = 1. Then there exists {D°'I°'v) {\x\k) = v{\x\k) for any x ^ 0. 

The proof consists of tedious but quite elementary calculations based on the integration 
formulas (2.2)-(2.8). A relatively nontrivial tool is the sum formula for the arithmetic-geometric 
progression (formula 0.113 from [Gradshteyn at al. 1996| ). ■ 

Using Lemma 3, we can consider the simplest Cauchy problem 



D'^uHxIk) = fi\x\ 
where / is a continuous function, such that 



m(0) = 0, 



X] Ifil^)] < if « 7^ 1, or l/l*?')! < if a = 1. 

l=m l=m 

The unique strong solution is u = I°'f; the uniqueness follows from the fact that the equality 
D'^u = (in the sense of V{K)) implies the equality u = const; see [Vladimirov et al. 1994] 
or [Kochubei 200 1]. Therefore on radial functions, the operators D" and I" behave like the 
Caputo-Dzhrbashyan fractional derivative and the Riemann-Liouville fractional integral of real 
analysis (see, for example, |Kilbas et al. 2006] ). However the next example illustrates a different 
behavior of the "fractional integral" in the non- Archimedean case. 

Example. Let /(|x|x) = 1, xeK. Then (I'^f) {\x\k) = 0. 

Proof Let \x\k = g". If a ^ 1, then by (3.10), (2.2), and (2.6), 

1 - g-" 



1 - 

\y\K<q 

1 -g-" 



\x 



K 



q""\x\K + 



a-l 



1 _ ^-1 
n-l| _ ^ ^a{n~l) 



1 - g-" 

g~"N^ + 



a Q 



-1 



If a = 1, then by (3.11), (2.4), and (2.6), 



(/V) (I 



x\k) 



g \x\k + 



glogg 



g mK + 



\y\K<q" 

n-1 



glogg 



(log|x|A- -log|2/|i^) dy 
1 



g" log \x\k — { n — 1 — 



g-1 



\x\k g 



-1 



i-g 

glogg 



1 



g-1 



g" ^logg 

g^Mogg 
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Of course, the above identities in the weaker sense of distributions from are trivial, 

since the constant functions are identified with zero, with and D°'D~°' — I. 

On the other hand, the example shows that the condition of decay at infinity in Lemma 3 
cannot be dropped. 



4 Fractional Differential Equations 

4.1. The Cauchy problem and an integral equation. In the class of radial functions 

u = u{\x\k), we consider the Cauchy problem 

D"u + a{\x\K)u^f{\x\K), xeK, (4.1) 

u{0) = 0, (4.2) 

where a and / are continuous functions, that is they have finite limits a(0) and /(O), as x — >■ 0. 

Looking for a solution of the form u — I°'v, where v is a radial function, we obtain formally 
an integral equation 

vHxIk) + a{\x\K) (rv) {\x\k) = f{\x\K): xEK. (4.3) 

Let us study its solvability. Later we investigate, in what sense a solution of (4.3) corresponds 
to a solution of the Cauchy problem (4.1)- (4.2). 

It follows from (4.3) that v{0) = /(O). Suppose first that a 7^ 1. By Lemma 2, the equation 
(4.3) can be written in the form 

[l + q-''a{\x\K)\x\'],]v{\x\K) 

+ cMk) J {\xr£'-\y\T')v{\y\K)dy^f{\x\K), x^O, (4.3') 



where c„ = 



\y\K<\x\K 

1-q-" 



1 - 

Since a is continuous, there exists such N & Z that 



q "'a{\x\K)\x\'^ < 1 for \x\k < q^ ■ 
On the ball Bn — \^x & K : \x\k < q^}, the equation takes the form 

v{\x\k)+ J ka{x,y)v{\y\K)dy = F{\x\k) (4.4) 

\y\K<\x\K 

where 

K{x,y)= [l + q-''a{\x\K)\x\%]-\M\x\K){\xrK-'-\y\%-'): x^Q, 

kaiO, y) = 0, F{\x\k) = [1 + q-''a{\x\K)\x\l\-^ /(kk). 
If we construct a solution of (4.4) on B^, and if 

a{\x\K) ^ for any X e X, m e Z, (4.5) 
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we will be able to construct a solution of (4.4), thus a solution of the equation (4.3), successively 
for all X & K. 

If q; = 1, we use (3.11) and obtain in a similar way the equation (4.4) with 

ki{x,y) = — [l + q'^a{\x\K)\x\K]~ a(|x|;^) (log - log |?/|x) , x 7^ 0, 

qlogq 

fci(0, y) = 0, F{\x\k) = [l + q-H\x\K)\x\Kr' /(kk). 

It is convenient to extend ka (including the case a = 1) by zero onto x B^. 

Theorem 1. For each a > 0, the integral equation (4-4) has a unique continuous solution on 

Proof. Let us consider the integral operator K, appearing in (4.4) as an operator on the 
Banach space C{Bj^) of complex- valued continuous functions on Bn- By the theory of integral 
operators developed in sufficient generality in [Edwards 1965] (see Proposition 9.5.17), to prove 
that /C is a compact operator, it suffices to check that, for any xq G -Bat, 

lim / \ka{x,y) - ka{xQ,y)\ dy = (4.6) 

X^Xq J 

Bn 

The relation (4.6) is obvious for xq 7^ 0, and also for a > 1. For xq = 0, we have A;q,(0, y) = 0, 
and for < a < 1, \x\k = q"', n> N, we get by (2.2) and (2.6) that 

\k4x,y)\ dy = const- J " 9"^""'^) dy 

= const ■ Q - = const -{xl^ 0, 

as \x\k 0. For a = 1, we use (2.4) and (2.6) to obtain that 

j \ki{x,y)\ dy = const- j (log g" - log |?/|if) d?/ 

Bm ls/|jf<g"-i 

= const ■ [nq^^^ log q — in — 1 )?"^^ log q 

. 9 n-i 1 log g I I „ 

= const q logg = const \x\k U, 

g — 1 — 1 

as \x\k — ^ 0. 

Thus, K, is compact, and by the Fredholm alternative ( [Edwards 1965] . 9.10.3), our theorem 
will be proved if we show that K, has no nonzero eigenvalues. 

Suppose that ICw = Xw, A 7^ 0, for some w G C{Bm). We have |w(?/)| < C, 

\k^{x,y)\<M\\x\l-^~\y\l7^^ 
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if a 7^ 1, and 

\ki{x,y)\ < M(log|x|i^ -log|?/|i^) 

if a = 1, \y\K < \x\k- 

In subsequent iterations we will deal with the integrals 



|j/k<k|jf 

and 

h,m= J i\og\x\K -\og\y\K)\y\Kdy. 

2/lA-<k'lK 

If a > 1, we find denoting \x\k = (t and using (2.2) that 



\y\T dy- J \y\T^'^-' dy = d^^MT^'^ 

y\K<q"-'^ \y\K<q"-^ 



where, for all m = 0, 1, 2, . . . 

' ]^ am— 1 ]^ q—am—a 



= (1 - 'Z"')?^ —TVr Z 7T < 

for some A > 0. A similar result, 

Ia,m = da^m\x\J \ da^rri ^ Aq , m = 0, 1, 2, . . . (4-7) 

is obtained for < a < 1, so that (4.7) holds for any a ^ 1. If « = 1, then the integral Ji m is 
evaluated as follows. We have 

n—l p ^ n— 1 

m+1) 



J2 J {^og\x\K-\og\y\K)mdy={l--)\ogq 

k=—cio, I I, ^ A:=— oo 

\y\K=q'' 

, oo 

= (1 - -) log g y = di 

q ^-^ 

where 



a(m+l) 



di,m = (1 - -) logqJ2 uq-"^""^'^ = (1 - -) logg ^ '"^ < Ag— 

g ^ g (i-g-"^-!)^ 

(we have used the identity 0.231.2 from [Gradshteyn at al. 1996| ). Thus, we have proved (4.7) 
also for a = 1. 

Let us return to a function w satisfying the relation K.W = Xw, A 7^ 0. Using (4.7) (separately 
for a 1 and a = 1) and iterating we find by induction that 

k(a;)| < C(M|ArM)'" ^fjg-"-'^ m = 0,1,2, x E B^. (4.8) 
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Since fj Q~"^ = g-fm(m+i)^ j^- follows from (4.8) that w{x) = 0. ■ 

j=0 

4.2. Strong solutions. Below we assume that the inequahty (4.5) is satisfied. Then, as 
we have mentioned, the solution v of (4.4) is automatically extended in a unique way from Bj^ 
onto K. The extended function v satisfies the equation (4.3). Therefore the function u — I'^v 
satisfies (4.1) in the sense of distributions from The initial condition (4.2) is satisfied 
automatically. 

Let us find additional conditions on a and /, under which this construction gives a strong 
solution of the Cauchy problem (4.1)-(4.2). Note that, by Lemma 3 and Theorem 1, a strong 
solution is unique in the class of functions u — I°'v where is a continuous radial function, 

oo 

such that ^ \v{q'')\ < oo for sum m G Z. 

l=m 

Theorem 2. Suppose that 

\a{\x\K)\<C\x\],"-', \f{\x\K)\<C\x\-^', £>0,C>0, (4.9) 
as \x\k > 1- Then u — I'^v is a strong solution of the Cauchy problem (4-l)-(4-2)- 

Proof. Let f (|x|a') be the solution of the equation (4.3') constructed above for all x & K 
(for X = 0, the integral in the right-hand side is assumed equal to zero). For \x\k < q'^ the 
existence of a solution v was obtained from the theory of compact operators; for larger values 
of \x\k we use successively the equation (4.3') itself. Denote 

Vm^ sup |^;(g"^)|. 

\x\K<q"' 

The sequence {Vm} is non-decreasing. 

As we assumed in Theorem 1 only the continuity of the coefficient a, we took N in such 
a way that the neighborhood = {x : \x\k ^ q^} was sufficiently small. Here we assume 
(4.5), so that we can take any fixed integer N and obtain a solution v on Bj^. 

Consider the case where a 7^ 1. It follows from (4.5) and (4.9) that 



[l + g-"a(kkkl^] ' 



< H 



where H > does not depend on x e K. li m > N, then we find from (4.3') and the above 
estimate for 1^^ that 

\v{qn\ < Cada,oHa{q"')q''"'V^., + H \f{g"^)\ . (4.10) 
Let us choose rui > N so big that 

c^d^,oHa{q^)q""^ <\, H \f{g^)\ < ^V^_^, 

as m > mi (that is possible due to (4.9)). Then it follows from (4.10) that Vm < Vm-i, as 
m > mi, hence that the function v is bounded on K. 
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Now we get from (4.3') and the assumptions (4.9) that 

\v{\x\k)\ < C\x\-^' , \x\k>1, (4.11) 

C > 0. A similar reasoning works for a = 1. 

Taking into account the estimate (4.11) we find from Lemma 3 that {D°'I°'v){\x\k) = 
v{\x\k)i X 7^ 0. Therefore the function u = I^v satisfies the equation (4.1) for all x 7^ 0. Since 
a, /, and u are continuous, the equation is satisfied in the strong sense. ■ 

4.3. Generalizations. Instead of (4.2), one can consider an inhomogeneous initial con- 
dition m(0) = Mo, Mo £ C. Looking for a solution in the form u = Uq + I°'v, v = v{\x\k), we 
obtain the integral equation 

v{\x\k) + a{\x\K) {I"v) {\x\k) = f{\x\K) - a{\x\K)uo, 

which can be studied under the same assumptions. 

All the above results carry over to the case of a matrix-valued coefficient a(|x|i^) and vector- 
valued solutions. In this case, to obtain a strong solution, it is sufficient to demand that the 
spectrum of each matrix a(|x|x), x E K, does not intersect the set {—q^,N G Z}. 
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